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Abstract 

The acoustic emission of fracture precursors, and the failure time of 
samples of heterogeneous materials (wood, fiberglass) are studied as a 
function of the load features and geometry. It is shown that in these 
materials the failure time is predicted with a good accuracy by a model 
of microcrack nucleation proposed by Pomeau. We find that the time 
interval St between events (precursors) and the energy e are power law 
distributed and that the exponents of these power laws depend on the 
load history and on the material. In contrast, the cumulated acoustic 
energy E presents a critical divergency near the breaking time r which is 
E ~ {■^-^) ■ The positive exponent 7 is independent, within error bars, 
on all the experimental parameters. 



1 Introduction 

Heterogeneous materials are widely studied not only for their large utility in 
applications but also because they could give more insight to our understand- 
ing of the role of macroscopic disorder on material properties. The statistical 
analysis of the failure of these materials is an actual and fundamental problem 
which has received a lot of attention both theoretically ^, ^, ||, 01 and 

experimentally!^, ||, |l^, ^ When an heterogeneous material is stretched 
its evolution toward breaking is characterized by the appearance of microcracks 
before their final break-up. Each microcrack produces an elastic wave which 
is detectable by a piezoelectric microphone. The microcraks constitute the so 
called precursors of fracture. It is very well known that these materials sub- 
jected to a constant load may break after a certain time, which is a function 
of the applied load. Many models have been proposed to predict this failure 
time, but the physical mechanisms remain unclear[^ ^, Very recently it 



has been propose d p^ , a model, which explains quite well the failure time 
of microcrystals [[ITJ and gels [|l^ submitted to a constant stress load. This 
model is based on the idea that a nucleation process of a microcrack has to 
take place inside the materials, in order to form the macroscopic crack. This 
nucleation process is controlled by an activation law, as the coalescence of a 
phase into another in a liquid-solid transition. Based on this prediction [|l6j, L. 
Pauchard et al. jlTj found that if a constant load is applied to a bidimensional 
microcrystal[p^ , it breaks after a time r given by the equation t = Toe^°/^ , 
where P is the applied pressure, and Tq and Pq are constants. Bonn et al. 
[ p!9| found a similar law for gels. Pomeau predicted that for three-dimensional 
microscopic systems the life-time should be: 

r = Toexp(^) (1) 



where Tq is a characteristic time and Pq a characteristic pressure, which mainly 
depend on the material characteristics, the experimental geometry and temper- 
ature. The idea, that the life time of a material might be due to a thermally 
activated process, has been proposed long time ago by Mogi |Q and Zhurkov 
[ p!3| . They got a different expression for r 

T = Toexp(^--^^ (2) 

This equation was accurately checked in many homogeneous materials and it 
shows a good agreement with experimental data jl^. However it has to be 
stressed that eq.0 and || are certainly first order approximation because they 
neglect the tensorial nature of crack perturbations and their long range interac- 
tions 1^. These ideas are quite interesting and it is important to check experi- 
mentally whether they can be applied in heterogeneous materials, such as fiber 
glass and wood pannels|Q. In two recent papers 0, we have shown that in 
these materials the microcracks, preceding the main crack form something like a 
coalescence around the final path of the main crack. The purpose of this paper 
is to investigate more deeply the behaviour of these materials, and specifically 
the statistical properties of fracture precursors and their relationships with the 
failure time. The paper is organized as follows: in section 2 we describe the 
different experimental settings. In section 3 we study the time failure for the 
samples submitted to a constant load (3.1) and its statistical properties (3.2). 
Then we generalize to a time dependent load. In section 4 the statistical be- 
haviour of the fracture precursors is studied as a function of the load features 
and the geometry. Discussion and conclusions follow in section 5. 

2 Experimental setups 



In order to verify the dependence of the results on the geometry and on the 
fracture's mode, we used three different experimental setups. We performed 
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mode-I fracture experiments both with a classical tensile machine (TM) and a 
high pressure chambers (HPC). In the case of the TM the stress distribution 
is very simple but, due to moving mechanical parts, we have to deal with a 
large acoustic noise. In order to avoid noise and to detect reliably microcracks 
with a weak sound emission, we have designed a set-up in which there are no 
moving parts except the sample itself, the HPC. Here, the stress distribution is 
very complicated but numerical calculations show that the experience can 
be thought of as a Mode-I test with circular symmetry. Finally, we used a 
flexion-machine to perform mode-I tests. 

2.1 Samples 

Several materials have been used. Most of the runs have been carried out on two 
fibrous composite materials: chipboard wood panel, which is made of small wood 
fibers randomly oriented, and different fiberglass panels made of a fiber fabric 
and an epoxy resin. The Young modulus of the samples is 1.8 10^ N/m? and 
10^° N/m^ for the wood panels and the fiberglass respectively. The Poisson's 
modulus isv^ 0.35 for both materials. The longitudinal sound velocity is 1900 
m/s for wood panels and 2200 m/s for fiberglass. The choice of the materials 
was determined by their features: they consist of small fibers, randomly oriented, 
and they are elastic and heterogeneous. The geometry of the samples, which 
depends on the experimental set up used to test it, is described in the following 
sections. 

2.2 High pressure chambers (HPC) 

A circular wood or fiberglass sample having a diameter of 22 cm and a thickness 
between 1 and 5 mm is placed between two chambers between which a pressure 
difference P = P2 — -Pi is imposed(see fig. 0b). If the deformation of the plate at 
the center is bigger than its thickness, which is the case here, the load is mainly 
radial |23[ |2^. Therefore, the experience can be thought of as a Mode-I test 
with circular symmetry. The pressure difference P supported by the sample is 
slowly increased and it is monitored by a differential transducer. This measure 
has a stability of 0.002 atm. The fracture pressure for the different tested 
materials ranges from 0.7 to 2 atm. We regulate P by means of a feedback loop 
and an electronically controlled valve which connects one of the two chambers 
to a pressurized air reservoir. The time taken to correct pressure variations 
(about 0.1 second) is smaller than the characteristic time of the strain rate. An 
inductive displacement sensor (Linear Differential Variable Transducer 500HR 
from PM Instrumentation) gives the deformation at the center of the plate 
with a precision of about 10 microns (the deformation just before fracture is 
of the order of one centimeter, depending on the material). The apparatus 
is placed inside a copper box covered with a thick foam layer to avoid both 
electrical and acoustic noise. Four wide-band piezoelectric microphones (Valpey- 
Fisher Pinducer VP-1093) are placed on the side of the sample (see fig. ^,b). 
The signal is amplified, low-pass filtered at 70 kHz , and sent to a digitizing 
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oscilloscope and to an electronic device which measures the acoustic energy 
detected by the microphones. The signal captured by the oscilloscope is sent 
to a computer where a program automatically detects the arrival time of the 
acoustic emissions (AE) at each microphone. If the signal is detected by more 
than two microphones, a calculation yields the position of the source inside the 
sample. A fraction of the detected events are rejected, either as a result of a large 
uncertainty of the location, or because they are regarded as noise. The mean 
standard error for the calculated positions is about 6 mm, which results mainly 
from the uncertainty of the arrival time. The electronic device that measures 
the energy performs the square of the AE amplitude and then integrates it over 
a time window of 30 ms, which is the maximum duration of one acoustic event. 
The output signal is proportional to the energy of the events, and its value is 
sent to the computer. The dynamic range for the energy measurement is four 
decades, and the device is adjusted in such a way that only the strong sound 
emitted by the final crack saturates it. The global results of the measurements 
are the following: a list of the positions of microcracks, the strain of the samples 
and the energy released as a function of the control parameter P. Further details 
of the setup are described elsewhere [|o[ |9| ^ p2|. 

2.3 Tensile Machine (TM) 

The experimental apparatus consists of a tensile machine (see fig. ^) which 
can apply a maximum force of about 23000 N. During the load we measure 
the applied force F, the strain, the AE produced by microcracks and the time 
at which the event was detected. The data acquisition set-up is the same used 
with the HPC. The samples have a rectangular shape of size I = 30 cm, w = 20 
cm and thickness of 2 mm. More details of the experimental setup can be found 

in m. 

2.4 Bending-machine (BM) 

The apparatus is a three-points flexion machine, flg.§. The rectangular sample 
lies horizontally with fixed edges and vertical load is imposed in its center. The 
size of the samples is / = 8 to 22 cm, w = 1 to 2 cm, and the thickness is 0.2 
cm. With this setup we only used fiberglass samples. We can load the sample 
up to 65Kg (that is the machine critical load) by minimum steps of 100 g. An 
inductive displacement sensor, similar to the one used in the HPC, has been 
used to measure the displacement of the center of the sample /. The sensor is 
connected to a computer that samples the signal at IHz. The failure time r 
is obtained by the analysis of the signal f{t); The uncertainty of r is then 0.5 
s. No acoustic emissions are measured with this apparatus. More details are 
described in ]2^ . 
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3 The failure time 



The aim of this section is to study the hfetime t of heterogeneous materials 
and to check whether eq.(l) could be useful in order to predict it. We first 
investigate r when samples are submitted to a constant load using the HPC 
and the BM. Then we consider the case of a time dependent load, and we try 
to generalize eq.(l). 



3.1 Constant load 

We first impose a constant strain to our samples (using all the apparatus), as 
it has been made for crystals As strain is fixed, every microcrack leads 
to a pressure decrease: in fact each microcrack weakens the material, so that 
a lower pressure is needed to keep the strain constant. It follows that if the 
imposed strain is small enough|^ the system reaches a stationary state, where 



the pressure remains constant and no more microcracks are detected 22 1 (see fig. 
||b). One sample was submitted to a large deformation (i.e. close to the critical 
value) and it did not break after three days. Therefore, at imposed strain, the 
effect observed in microcrystals is not valid for heterogeneous materials. 

On the other hand, if a constant stress is imposed to the system, no matter 
which apparatus we use, it will break after a certain time which depends on the 
value of the applied load. This can be done by imposing a constant pressure 
with the HPC (see fig. ||a) or a constant force with the TM and BM. The reason 
for this is that after every single microcrack the same load must be endured by 
the weakened sample, so that it becomes more and more unstable. Using either 
the TM or HPC, we have submitted several samples to different constant loads 
P and we have measured the life-time t. The values obtained are well fitted 
by eq.(|l|), that is the exponential function predicted by Pomeau. On the other 
hand, the life-time expression r = ae~^^ proposed by Mogipfl does not conform 



to our data|27 . The same law has been found experimentally by Zhurkov]13|. 
However it's worth noting that his work deals mainly with homogeneous, visco- 
elastic and plastic materials, whereas the materials we used are heterogeneous 
and elastic: this could explains why he found a different dependence of r on the 
imposed load. In fig. t is plotted versus -p^ in a semilog scale and a straight 
line is obtained. Fig .p[ corresponds to the case of wood samples broken in the 
HPC, while in fig. Hwe show the points for fiberglass samples broken in the 
BM. Each point corresponds to the mean value obtained with 20 samples. Even 
if the load is very small the sample will eventually break, although the life-time 
can be extremely long. For example, using eq.(|l|) and the best fit parameters 
of fig. one estimates r ~ 5000 s at P = 0.43 atm. Halving the imposed 
pressure causes r to become extremely large : t = 4.4 • 10'^^ years at P = 0.21 
atm). 

The value of Tq seems to depend on the geometry and the material but not, 
within the error bar, on the sample size. In fact, for the circular samples broken 

^it's clear that if the imposed strain has to be smaller than a critical value, at which the 
sample breaks instantaneously. 
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with the HPC, we find Tq = 50.5 ± 0.2 s for wood and To = 44.6 ± 0.2 s for 
fiberglass. While in the the BM we find To = 2.5 ± 0.3 s and = 2.7 ± 0.3 s 
for the samples with W = 1cm and W — 2cm respectively. 

The value of Po depends both on the geometry, the size and the material of 
the sample, and indeed we find Po = 2.91 atm and Po = 0.63 atm for fiberglass 
and wood broken in the HPC. For the fiber glass samples broken with the BM 
the values of Po are Po = 71.1 kg for W = 2 cm and Po = 35 kg ior W = 1 cm 



3.2 Statistic of the failure times 

It's interesting to study the statistical distribution N{t) of the failure times. 
This information can give more insight on the physical phenomena. 

The main hypothesis of Pomeau is that the failure of a sample is due to 
the thermal nucleation of one defect. Thus one expects that the failure time 
T follows a Poisson's distribution. This has been experimentally observed in 
gels[^. Conversely, in all the measures made with the HPC and the TM the 
failure time follows a normal distribution. This is also the case for crystals [|l7| . 

This difference can be explained by the fact that in our samples the failure is 
due to the nucleation and coalescence of a large number of defects, each of one 
is thermally activated and would eventually follows a Poisson's law, if it were 
isolated. Numerical simulations and analytical calculations seem to confirm this 
ideafU, |2|, |3^. 

In the case of the experiments performed with the BM, the failure time dis- 
tribution N{t) seems not to follow a Poisson's law nor a normal distribution, fig. 
In this case we observed that the cumulative distribution 

is best fitted by the sum of two exponentials, fig.Pa. We believe that this is 
due to the fact that the two components of these samples (the resine and the 
texture of glass fibers) give rise to different characteristic times. We think that 
this "separation effect" is observed only with the BM because of the small size 
of the samples. 



3.3 Time dependent load 

In order to find a law that holds for a time dependent imposed stress, we intend 
to generalize the eq.(0). If the pressure P changes with time, it is reasonable to 
consider the entire history of the load. Therefore we consider that 



— exp 

To 



is the density of damage per unit time, where Tq and Po are fitting parameters 
obtained in the constant load case. The certitude of breaking is obtained after 



6 



a time r such that: 

dt = 1 (3) 

where Tq and Po have the previously determined value. Notice that this equation 
is equivalent to eq.(^ when a constant pressure is applied. 

To check this idea, we have applied the load to the sample (using the HPC) 
following different schemes. We have first applied successive pressure plateaux 
in order to check whether memory effects exist. In fig. |^a the pressure applied 
to the sample is shown as a function of time. A constant load has been applied 
during a certain time ti, then the load is suppressed and the same constant 
load is applied again for a time interval T2 . The sample breaks after a loading a 
time Ti + T2 which is equal to the time needed if the same load had been applied 
continuously without the absence of load during a certain interval. Therefore a 
memory of the load history exists. The life-time formula (eq. ||) is also valid 
if different constant loads are applied successively (fig. |^b). This concept can 
explain the violation of the Kaiser effect in these materials [pO|. 

If the load is not constant, the life-times resulting from the proposed integral 
equation are still in good agreement with experimental data. A load linearly 
increasing at different rates Ap has been applied to different samples. The 
measured breaking times are plotted in fig. |^ along with a curve showing the 
values computed from eq.^ Even if a quasi-static load is applied erratically (fig. 

, the calculated life-time agrees with the measured one. These experiments 
show that eq. ^ describes well the life-time of the samples submitted to a time 
dependent pressure. 




3.4 The dependence of r on the temperature 

The question is to understand why eq.(|l|) and (^) works so well for a three 
dimensional heterogeneous material. Indeed, in the Pomeau formulation [ p^ 

^>.^g(!L_) (4) 

where Y is the Young modulus, T the temperature, K the Boltzmann constant 
and T] the surface energy of the material under study. G is a geometrical factor 
which may depend on the experimental geometry, on defect shape and density. 

In our experiment with HFC, we found Po = 0.62 atm for wood, which has 
Y=1.8- 10® TV/m^, and Po = 2.91 atm for fiberglass, which has Y = 10^° N/m^. 
Thus the ratio between the values of Po found for the two materials is closed to 
the ratio of the square root of their Young modula. 

In contrast temperature does not seem to have a strong influence on r. 
In fact we changed temperature, from 300K to 38Gi^ which is a temperature 
range where the other parameters, Y and rj, do not change too much. For this 
temperature jump one would expect a change in r of of about 50% for the 
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smallest pressure and of about 100% for the largest pressure. Looking at fig. || 
we do not notice any change of t within experimental errors which are about 
10%. In order to maintain the change of r within 10% for a temperature jump 
of 80K one has to assume that the effective temperature of the system is about 
3000K. Notice that this claim is independent on the exact value of the other 
parameters and G. 

These observations seem to indicate that the nucleation process of microc- 
racks is activated by a noise much larger than the thermal one. Such a large 
noise can be probably produced by the internal random distribution of the de- 
fects in the heterogeneous materials that we used in our experiments. This 
internal random distribution of material defects evolves in time because of the 
appearance of new microcracks and the deformation of the sample. Therefore 
this internal and time dependent disorder of the material could actually be the 
mechanism that activates the microcrack coalescence and play the role of a very 
high temperature. Numerical simulations and analytical calculations, which we 
performed in fuse networks, confirm this hypothesis p8|, p9| p0[ . Similar conclu- 
sions on the role of disorder in the crack activation processes have been reached 
by other authors]^, The experimental test of these models is a very impor- 
tant point which merits to be deeply explored in the future. 



4 Statistical behavior of fracture precursors 

When a constant pressure is applied to the sample using the HPC, the acoustic 
emissions of the material are measured as a function of time. We find that the 
cumulative acoustic energy E diverges as a function of the reduced time 
specifically E oc {^^y with 7 = 0.27 (see Fig. |l^). Notably, the exponent 7, 
found in this experiment with a constant applied pressure, is the same of the 
one corresponding to the case of constant stress rate [^. Indeed it has been 
shown [|[ 10) that if a quasi-static constant pressure rate is imposed, that is 
P = Apt, the sample breaks at a critical pressure Pc and E realesed by the final 
crack precursors (microcracks) scales with the reduced pressure or time (time 
and pressure are proportional) in the following way: 

where r = Pc/Ap in this case. Thus it seems that the real control parameter of 
the failure process is time, regardless of the fact that either a constant pressure 
rate or a constant pressure is applied. In the case of constant load rate {P = Apt 
OT u = Bt) the system has not a characteristic scale of energy or time: the 
histogram N{e) of the released energy and the histogram N{St) of the elapsed 
time St between two consecutive events reveal power laws, i.e. N{e) ~ and 
N{5t) ~ St~^ . The exponents a, jS and 7 do not depend on the load rate Ap or 
|l^. Power laws for similar magnitudes are found experimentally on cellular 
glass |l^ , and numerically in a related process, the dielectic breakdown p2[. The 
value of the exponents are not too different. We are interested in studying the 
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exponents in different geometries and wlien a constant (creep test), cyciic or 
erratic load are imposed. To check the dependence of a, (3 and 7 on the geometry 
we used the TM. The force apphed to the sample is slowly and constantly 
increased till the sample fails. During the loading we measure the applied force 
F , the strain, the AE produced by microcracks and the time at which the event 
was detected. 

4.1 The dependence of a and (3 on the load features 

In the experiments performed with the HPC, power laws are obtained for the 
distributions of e and of 5t. As an example of two typical distributions obtained 
at constant imposed pressure, we plot in fig ^) and ^) N{5t) and N{5e) re- 
spectively. The exponents of these power laws (etc for energies and (3c for times) 
depend on P. In fig. |^, ac and (3c are plotted versus P. Note that both expo- 
nents grow with pressure. We observe that the rate of emissions increases with 
pressure, so that the weight of big values of 5t decreases. This explains the fact 
that [3c grows with pressure. We have compared the histograms of energy e for 
several pressures, and we noticed that the number of high-energy emissions is al- 
most the same, while the number of low-energy emissions increase with pressure, 
so that the exponent ac increases as well. Moreover, as the pressure increases, 
the exponents ac and (3c attain the values a = 1.9 ± 0.1 and /? — 1.51 ± 0.05 
obtained in the case of a constant loading rate[|l^. We imposed to the sample 
a erratic and an cyclic load, which are plotted as a function of time in figure 
prj a and p] b respectively. Power laws are obtained for the distributions of e and 
for 5t. The exponents of these power laws do not depend on the load behavior; 
their value is the same of that at constant loading rate. These and previous 
results H, |o|, allows us to state that if ^ 7^ 0, the histograms of the released 
energy e and of the time intervals 5t do not depend on the load history. The 
fact that a and (3 do not depend on ^ seems to be in contrast with the fact 
that ac and (3c depend on P. This result can be interpreted by considering that 
the microcracks formation process is not the same when ^ — and ^ 7^ 0. In 
the former case, imposed constant P, the mechanism of microcrack nucleation 
is the dominant one and the nucleation time depends on pressure. In the other 
case, ^ 7^ Oi the dominant mechanism is not the nucleation but the fact that, 
when pressure increases as a function of time, several parts of the sample may 
have to support a pressure larger than the local critical stress to break bonds. 
The fact that at high constant pressure ac and Pc recover the value ac and (3c 
has a simple explanation. Indeed, in order to reach a very high pressure Ph, ^ 
is different from zero for a time interval which is comparable or even larger than 
the time interval spent at constant pressure P^. Thus at high constant pressure 
the system is close to the case ^ 0. 

4.2 The dependence of 7 on the load features 

The measures performed with the HFC imposing an erratic and a cyclic pres- 
sure, plotted respectively in fig. pT| a and pT[b, allow us to check the dependence 
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of 7 on the history of the sample, i.e. on the behavior of the imposed pressure. 
The cumulated energy E for the erratic and the cyclic pressure, shown in fig 
|n|a and ^ as a function of t, is plotted in log-log scale as a function of the 
reduced parameter ^^-^ in fig and [l^ b respectively. We observe that, in 
spite of the fluctuations due to the strong oscillations of the applied pressure, 
near the failure the energy E, as a function of 2-^, is fitted by a power law 
with 7 ~ 0.27 ± 0.02. In fig.p^ (reproduced in fig.^ for the sake of clearness), 
E measured when a constant pressure is applied to the sample is plotted as a 
function of -^-y^- A power law is found in this case too The exponent 7 is, 
within error bars, the same in the three cases. Hence it seems not to depend 
neither on the applied pressure history nor on the material lol 



Further, experiments made with the TM show that 7 is independent on the 
geometry. In fact we observe that the behavior of the energy near the fracture 
as a function of (-^) is still a power law of exponent 7 ~ 0.27, as shown in 
figure Hd. 



5 Discussion and conclusions 

We presented the results of experiments regarding fracture of heterogeneous 
materials (chipboard and fiberglass). When these materials are submitted to 
a strain, one observes acoustic emissions (precursors) before the samples fails. 
We have measured the acoustic emissions and the lifetime t of the samples. 

We have shown that eq.(l) proposed by Pomeau predicts well the functional 
dependence of r on the applied load. However the original Pomeau's theory 
is unable to explain some features that are observed in our experiments, in 
microcrystals [p7| and gels|jl^. In fact eq.(l) is based on the idea that the 
fracture is due to the nucleation of one preexisting defect, which is thermally 
activated. We have shown that in our experiments the fracture is due to the 
nucleation and coalescence of a large number of defects, as confirmed by the 
presence of acoustic emissions and the shape of the statistical distribution of 
lifetimes. Moreover, we found that the lifetimes are not affected, within the 
experimental errors (20 %), by the temperature. We have calculated that in 
order to estimate the measured lifetimes, the temperature T to insert in eq.|] 
and in eqM should be about 3000 K. Similar results have been found in 2D 
crystals ]18[ and gels, where the temperature T to insert in eq.(3) should be 
about lOOO/C < T < 2500ii: and T^// > 10^° K respectively. As for wood and 
fiberglass, also in gels the lifetime r of the sample is not influenced, in the limit 
of experimental errors, by a variation of the temperature T from 20 to 90 °C. 
In contrast, experiments on 2D-crystals show that r depends on T. 

To explain these points we propose to modify in a statistical model the 
one of Pomeau. In this revised model the failure of the sample is due to the 
nucleation and coalescence of a certain number of defects, each of which is 
thermally activated. The parameters Y, -q and T of eq.(3), become average 
parameters, which keep into account that the interaction between defects have 
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a tensor nature and are long range To explain the fact that the thermal 
temperature has a minor influence on the lifetime, we suppose that the strong 
time-dependent fluctuations of the internal forces induced by the heterogeneity 
(defects, microcraks ...) can be considered as a sort of noise. Thus T depends 
on the thermal temperature but mainly on the disorder in the medium. In 
this way the heterogeneity of the material enhances thermal fluctuations so 
that the nucleation time of defects becomes of the order of the measured ones. 
Recently Arndt et al. reached a similar conclusion [^if using a different approach, 
which is coherent with the generation of a time dependent distribution of the 
microcrack ensemble . Numerical simulation that we have performed on a very 
simple model|^, ^ are in agreement with these results. This model allows us 
to generalize eq.(l) to the case of a time dependent load, and to explain the 
violation of the Kaiser effect in these materials. 

We also studied the statistical properties of the fracture precursors. We have 
found that the histograms of the energy and of time between two consecutive 
AE follow power laws of exponents (3 and a respectively. In proximity of the 
fracture the cumulated energy follows a power law, typical of phase transitions. 
Notably, the critical exponent 7 seems to be independent on the geometry and 
the applied load. Indeed if AE is considered as a susceptibility it is not easy to 
put together the observed critical divergency with a nucleation process. Prob- 
ably the standard phase transition description can be only partially applied to 
failure because of the intrinsic irreversibility of the crack formation. 
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Figure 1: a, b) Sketch of the high pressure chamber (HPC) apparatus. S is the 
sample, DS is the inductive displacement sensor (which has a sensitivity of the 
order of 1 ^m). M are the four wide-band piezoelectric microphones. P=Pi-P2 
is the pressure supported by the sample. P is measured by a differential pressure 
sensor DPD( sensitivity 0.002 atm). EV is the electronic valve which controls 
P via the feedback control system Ctrl . c) Sketch of the tensile machine. 
An uniaxial force, which is measured by a piezoresistive sensor, is applied to 
the sample by a stepping motor. Four wide-band piezoelectric microphones 
measure the acoustic emissions emitted by the sample. Experiments have been 
done using rectangular (20 x 29 cm) wood samples of 4 mm thickness. The 
whole apparatus is surrounded by a Faraday screen. 
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Figure 2: The bending machine (BM). A vertical force is imposed in the center 
of the sample by using weights (up to 65 Kg)^ so that the sample is broken by 
bending. The edges of the sample can be clamped or free. The displacement is 
measured by the sensor DS. 
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Figure 3: A sample is submitted to a constant load with the HPC. During 
the load we measure the pressure (continuous line) and the deformation of the 
sample in its center (dashed line), (a) A constant pressure is imposed to the 
sample. The deformation increases continuously - even after that the pressure 
has reached a constant value - till the sample fails, (b) A constant deformation 
is imposed to the sample. In this case, after a transient period, the pressure 
decreases till the system reaches a stationary state (not shown in this picture). 
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Figure 4: Measurements on wood samples. The time r needed to break the wood 
samples under an imposed constant pressure P is here plotted as a function of 
1/P'' in a semilog scale. The dashed line represents the solution proposed by 
Mogi [^(t = ae^^^). The continuous line is the solution proposed by Pomeau 
for microcrystals (r = Tq^'^^"^^^ ). In the plot Tq = 50.5 s and Pq = 0.63 
atm. Every point is the average of 10 samples. The error bar is the statistical 
uncertainty. For the fiberglass samples, we find Tq — 44.6 s and Po — 2.91 atm. 
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Figure 5: Failure time r of the samples in fiberglass broken with the FM 
(clamped edges). The sample's size are 22 x 2 x 0.2 cm (a) and 22 x 1 x 0.2 cm 
(b). Each point represents the mean value of 20 measures. Lines represent the 
best fit with t — Tq exp (-^) (solid line), r = Tq cxp (^) (dashed-dotted line), 
r = ylexp(-6F) (dotted Hne), and r = A P^^ (dashed line). 
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T (S) 

Figure 6: Distribution of the lifetimes r of 100 samples in fiberglass broken 
with the BM (clamped edges) with a load P — 54 Kg. The sample's size are 
22 X 2 X 0.2 cm. a) The histogram of t shows that the distribution of lifetimes 

is not gaussian. b) The cumulative distribution Q{t) = jgo^^^*^^^*^ (solid line) is 

best fitted by the sum of two exponential terms (dotted line). 
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Figure 7: The imposed time dependent pressure (bold dotted line) is plotted 
as a function of time in the case of wood samples. The continuous line is the 
integral in time of the function / (P) = ■^e~^°/-^'^ . On the basis of eq. || the 
predicted breaking time r is obtained when the integral of / (P) is equal to 
1. The horizontal distance between the two vertical dashed lines in each plot 
represent the difference between the predicted and the measured breaking time. 
In a) a constant pressure has been applied during about 700 s, then the load 
is suppressed and then the same constant load is applied again. The difference 
between the life-time predicted by (eq. ||) and the experimental result is of 3%. 
b) Here two pressure plateaux of different value are successively applied to the 
sample. The difference between the measured and the predicted life-time is of 
5%. In c) an erratic pressure is applied to the sample. Here the error is of 10%. 
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Figure 8: A load linearly increasing at different rates Ap has been applied to 
different samples. The measured breaking times are plotted as a function of 
Ap in a loglog scale; circles and squares represent the measures on wood and 
fiberglass samples respectively at T=300 K. Bold circles represent measures on 
wood samples at T=380 K. The lines are the life time calculated from eqj| using 
the best fit values for Po and Tq. These experiments show that eq. ^ describes 
well the life-time of the samples submitted to a time dependent pressure. 
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Figure 9: a,b) Two typical time St and energy e distributions obtained at 
imposed constant pressure {P = 0.56 atm). c) The exponents ac (empty circles) 
and (3c (black points), plotted as a function of the value of the imposed constant 
pressure. Note that as the pressure increases, the values of the exponents tend 
to those obtained in the case of constant pressure rate. The error bars represent 
the statistical uncertainty. 
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Figure 10: The cumulative energy E, normalized to Emax, as a function of the 
reduced control parameter at the neighborhood of the fracture point (Case 
of imposed constant pressure). The circles are the average for 9 wood samples. 

The soHd hne is the Ht E = Eq (^-7^) ^. The exponent found, 7 = 0.26, does 
not depend on the value of the imposed pressure. In the case of a constant 
pressure rate the same law has been found 110, M . 
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Figure 11: The imposed pressure, normalized at Pi„ax-(^olid line) and the 
cumulative energy E, normalized to Emax, (dashed line) are plotted as a function 
of time fa) An example of erratic pressure, b) A cyclic pressure. 



24 



(x-t)/T 



(T-t)/T 




0.01 10 1.00 0.01 0.10 1.00 

(T-t)/T (^-t)/-C 



Figure 12: The cumulative energy E, normalized to ii^max, as a function of the 
reduced control parameter at the neighborhood of the fracture point. Figure 
d) represent the measure taken, at imposed constant rate force, on the tensile 
machine. The other figures represent measures made on the HFC apparatus at : 
imposed constant pressure (c), imposed cyclic pressure (a) and imposed erratic 
pressure (b). The dotted lines are the fit E = Eq (^-^) The exponents found 
are: 7 = 0.29 (a), 7 = 0.25 (b), 7 = 0.29 (c) and 7 = 0.27 (d). In the case of 
a constant pressure rate (on the HPC machine) the same law has been found 
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